Mixing and condensation in a wet granular medium 
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We have studied the effect of small amounts of added liquid on the dynamic behavior of a granular 
system consisting of a mixture of glass beads of two different sizes. Segregation of the large beads 
to the top of the sample is found to depend in a nontrivial way on the liquid content. A transition 
to viscoplastic behavior occurs at a critical liquid content, which depends upon the bead size. We 
show that this transition can be interpreted as a condensation due to the hysteretic liquid bridge 
forces connecting the beads, and provide the corresponding phase diagram. 

PACS numbers: 05.65.+b,45.70.-n,45.70.Mg 



As is generally known from everyday experience, the 
mechanical properties of a granular material change dra- 
matically if some liquid is added. The main reason is 
the internal cohesion due to capillaryforces arising from 
liquid bridges between the grains 
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recent years have seen considerable progress in under- 
standing the dynamics of dry granular materials U 0, 0| j 
the physical mechanisms underlying the properties of wet 
systems remain largely obscure. This is in part due to the 
inherent complexity of these highly nonlinear systems, 
but also to the fact that only few experiments have yet 
been conducted on this topic. 

Some experiments have been recently done on the in- 
fluence of liquid on the mixing behavior of systems con- 
sisting of large and small spheres [T(l ITll |. It was ob- 
served that as liquid was added, the tendency to segre- 
gate, which is well known for dry granulates 12, .l^l. was 
diminished. At a certain wetness, complete mixing oc- 
curred. In these experiments the morphology of heaps 
formed by pouring mixtures of particles with two differ- 
ent sizes through an orifice was investigated. In order to 
have a closer look at the observed effects, we have carried 
out experiments in which the agitation of the material is 
provided in a more homogeneous way, thus ensuring a 
more or less well defined granular temperature. Further- 
more, the propensity of the system to segregate into large 
and small beads can be determined quantitatively in our 
setup. We observe a strong dependence of mixing on the 
liquid content, and the transition to a viscoplastic regime 
due to condensation of the grains into clusters, which was 
reported in ref.j^, was clearly observed. However, the 
overall behavior was found to be more complicated than 
previously reported, and we are able to explain all of 
the observed features within a rather simple theoretical 
model. 

A typical sample consisted of 150to^ of glass beads with 
R = 2.bmm radius, sealed in a cylindrical jar (diameter: 
80 mm, cf. fig. 1) together with the same volume of 
smaller glass beads with radius r, which was varied from 
50 to bOOfim, and a well defined amount of liquid. In or- 



der to prevent crystallization, the small beads were cho- 
sen to be slightly polydisperse: the spread in bead size 
ranged from 10 to 20 %. The jars were mounted on a 
shaker, with the cylinder axis upright. The shaking mo- 
tion consisted of small horizontal circles with a diameter 
of 5 mm and a frequency of 20 revolutions per second, 
which corresponds to a velocity of the container walls of 
0.38 m/s. When the segregation pattern did not change 
any more (typically after one half hour), photographs of 
the jars were taken and analyzed. Two examples are 
shown in Figure 1. 




'Electronic address: Dimit rios. GeromichalosQphysik.uni- ulm. de| 



FIG. 1: Horizontally shaken beads (with radii R — 2.5mm 
and r = 0.5mm) without water (left) and after adding 75 /il 
(right). 

We defined the mixture zone as the zone in which both 
the smaller and the bigger particles were visibly present. 
This zone was observed to be rather homogeneous con- 
cerning the mixing ratio. Occasional analysis of horizon- 
tal slices of the samples revealed a small radial variation 
of the mixing ratio, which is not captured by the pho- 
tographs, but the results given below are qualitatively 
valid for the whole volume. For the results presented 
here we used water as the liquid, but similar results were 
obtained using ethanol. We can therefore exclude that 
leaching of ions from the glass, and other side effects due 
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to the specific properties of water, play a significant role 
[l^ . The quantity we will discuss is the height of the mix- 
ture zone, divided by the total height of the sample, and 
is denoted by h. For perfect segregation, h ^ 0, while 
h = 1 for perfect mixing. The amount of liquid added to 
the samples is denoted by the dimcnsionless quantity W, 
which is defined as the volume of the liquid divided by 
the total volume of the small beads (including the space 
between the beads). 

In contrast to what one might anticipate on the basis 
of earlier work 0, , the propensity of mixing of the 
particles did not improve monotonically with the amount 
of liquid added. Instead, we found in most cases that the 
height of the mixture zone, h, was first reduced upon 
addition of water. Only after a certain amount of liquid 
had been added did the mixing improve. In Fig. 2, h is 
plotted as a function of W for two different r. The size 
of the large beads was kept constant. Clearly, there is a 
rather well defined liquid content above which the height 
of the mixing zone increases when more liquid is added. 

The boundary between these two regimes (i.e., the 
minima of the curves as in fig. 2) is indicated by the 
open squares in Figure 3, which represents a phase dia- 
gram of the system in the plane spanned by the radius 
of the small beads, r, and the liquid content, W. The 
dotted line serves as a guide to the eye. The full circles 
correspond to the occurrence of complete mixing (i.e., 
h — I). We can thus clearly distinguish three regimes, 
which we denote as gaseous, intermediate, and viscoplas- 
tic 15, 16], respectively. In what follows, we will discuss 
these regimes in some detail. 



Segregation phase diagram 
(Big bead radius: 2.5 mm) 





0.006 



0.000 0.001 0.002 0.003 0.004 0.005 0.006 
W 



FIG. 3: Segregation phase diagram in dependence of the liq- 
uid added. In the grey area there occurs no segregation, while 
in the white area segregation was observed. The dashed steep 
line separates the area in which the segregation increases by 
adding water (left) from the area in which the segregation 
decreases or stays the same (right). 



tion, which occurs as ft. ^ 1, appears quite abrupt when 
inspecting the experiment by eye. While the general mo- 
tion of the whole sample was a kind of circular convection 
around the axis of the jar in the gaseous and intermedi- 
ate regimes, the motion in the case of perfect mixing 
was reminiscent of kneading raisin dough. We identify 
this transition with the one reported recently by other 
authors 0, We assume that this transition into 

the viscoplastic regime occurs when the granular tem- 
perature of the system, i.e., the average kinetic energy 
per bead, is not anymore sufficient to break the liquid 
bridges between the beads, and is thus a kind of conden- 
sation within the granular gas, mediated by the liquid 
bridges. 

In order to support this view, we have to discuss the 
involved forces in some detail. The force exerted upon 
two adjacent beads by the liquid bridge can be written 
as F{d,w) — 2nrjf{w)g{d), where 7 is the surface ten- 
sion of the liquid, r is the radius of the spheres, w is the 
liquid volume of the bridge, and / is some function de- 
pending on the roughness of the sphere surfaces 0, ^3 ■ 
g{d) describes the dependence of the force on the mutual 
distance, d, of two ideally spherical beads. The energy 
required to destroy the bridge is given by 



FIG. 2: Normalized height of the mixture zone of shaken 
glass beads as a function of the amount of water added. The 
radius of the big beads was R — 5mm, the diameter of the 
small beads was r = 240/im (closed symbols) and r = 140/im 
(open symbols), respectively. 

Let us start with the transition to the viscoplastic 
regime, in which perfect mixing is observed. This transi- 



Eb{w) = 27rr7/(w) / g{d)dd 



(1) 



where dp is the distance at which the bridge pinches off. 
The integral can be well approximated as [T^ . 

For perfect spheres and complete wetting, f{w) = 
1 Vw, but for a rough surface, it tends to zero as w — > 0. 
The shape of f{w) has been discussed in detail before 
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0. Based on these considerations, it can be well approx- 
imated by 



with Wo := rb"^ characterizing the roughness amplitude, 
8. For the exponent, we have fx = {2 — x)/(2 + x), where 
X is the roughness exponent of the bead surface. For 
(typical) scratch and dig roughness, x is close to zero, 
and n ~ I. This corresponds to the regime found by 
Hornbaker et al. 3]. In the opposite case, when x — > 1, 
we have — > 1/3. We see that /i is confined to a narrow 
interval, and is very probably close to one in our system. 

We now must calculate the average bridge volume, w, 
from the dimensionless liquid volume, W. At equilib- 
rium, assuming that all of the liquid is in the bridges, 
w = IlllJil where c is the coordination number of the 
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network of bridges (c ~ 6) and pp is the packing density 
of the spheres {pp w 0.65). Under dynamic conditions the 
typical bridge volume would be expected to be somewhat 
smaller than the equilibrium value. Our measurements 
[2^ suggest that w = ar^W, where a « 0.25. If we now 
equate the bridge energy, Eb, to the average kinetic en- 
ergy of the small beads, Ek ■= ^ (^^)' '^^ obtain the 
implicit equation 

^2 _ ;1/p^1/3m^1/3m^(2-1/p) + ^ Q (3) 

aW 

where we have introduced the characteristic length scale 
I := ^^.jj^ (p is the mass density of the bead material. In 

our case, p — 2500kg/m'^). Eq. (3) can be easily solved 
analytically for p. = 1 and p — ^. For the most relevant 
case of /Lt = 1, we obtain 

^-r'"W"' Jl + y'l-(^)"'J (4) 

which intersects the W-axis at the critical moisture con- 
tent Wc = ^(2(5/Z)^/^, which only depends on the rough- 
ness amplitude, (5, but not on r. This function represents 
the phase boundary at which condensation sets in, and is 
indicated as the full curve in Figure 3. From this best-fit 
curve we obtain S ~ 500nTO, which is similar to the peak- 
to-peak roughness obtained from the inspection of the 
beads by atomic force microscopy. From the fitted value 
of the intrinsic length scale I , and from the known density 
of the glass and the surface tension of the w ater, we get 
the average granular velocity, v := \/{v^ = O.llm/s. 
This compares favorably with the velocity of the shak- 
ing motion, which is expected to be of the same order of 
magnitude. 

Let us now turn to the gaseous regime. If there is no liq- 
uid in the sample at all, segregation is observed as usual, 
which has become known as the 'brazil nut effect' (BNE) 
El El 113 113. As it has been shown recently |3, this 
may be understood in a continuum picture, as follows. A 



granular medium is strongly dissipative due to the inelas- 
tic character of the impact kinetics of the grains, which 
is characterized by the restitution coefficient, e. The lat- 
ter is defined as the ratio of the particle velocities after 
and before impact, e = v'^^Jvrei- Thereby the granular 
temperature bath, which typically corresponds to ther- 
mal temperatures in the Terakelvin range, is internally 
coupled to the room-temperature bath provided by the 
atoms constituting the grains. This inherent state of non- 
equilibrium is the basic reason for the peculiar features 
of granular matter dynamics in general. 

Now consider a spherical region with radius R within 
the granular gas. When a grain, or small bead, enters this 
region, it will have many impacts, and thereby lose a lot 
of energy, before it can leave this region again. When 
the latter is replaced by a solid sphere (i.e., a large bead 
of radius R), the grain will experience only one such im- 
pact. Consequently, the granular gas is hotter close to an 
intruder than far away in the bulk, because the intruder 
is a less effective heat sink. The granular gas is therefore 
also less dense around the intruder, and the latter rises 
due to the buoyancy of the 'bubble' of reduced density it 
carries around itself E3- 

The height of the mixing zone is then a direct conse- 
quence of the balance between the hydrostatic pressure 
exerted by the large-beads fraction upon its lower inter- 
face, and the buoyancy pressure exerted upon the same 
interface from below by the mixed phase. This picture 
is supported by the observation that the large beads ac- 
quire a rather uniform spacing within the mixed phase, 
which might then correspond to (twice) the thickness of 
the reduced density layer around each large bead [l3j |. 

The BNE increases as the energy loss upon impact is 
increased. This, however, is just what is achieved when 
liquid is added. The interparticle force exerted by the liq- 
uid has a hysteretic character |17il22i.l2^ . psf and is thus 
intrinsically dissipative. We thus expect the restitution 
coefficient to be reduced, and the BNE to become more 
prominent, when liquid is added. More precisely, the 
buoyancy of the reduced density bubbles around the large 
beads increases, and a shallower mixing zone is needed to 
balance the hydrostatic pressure of the large-beads col- 
umn above it. Thus the height of the mixing zone is 
expected to decrease when liquid is added, which is what 
we observe. We thus assume that throughout this regime, 
which is defined by a monotonically decreasing height of 
the mixing zone, the physics of two-particle impact is af- 
fected by the liquid merely in the sense of a decrease in 
the restitution coefficient, e. 

Although a quantitative prediction of the BNE is yet 
difficult, we can give a rough estimate of the expected 
size of the effect. The inelastic character of the impact of 
the dry particles reflects itself in the BNE (and thus the 
tendency of demixing) of the dry material. Consequently, 
we expect the contribution due to the liquid to be of the 
same order of magnitude when the energy stored in the 
liquid bridges is similar to the energy lost upon impact, 
i.e., when Eb « (1 — e^)EK- Since I oc E]^^ , the corre- 



4 



spending moisture content, Wbne, is related to Wc via 
Wbne « (1 - e^)^/^tVc- Exploiting that Wbne roughly 
corresponds to the open symbols in Fig. 3, we can esti- 
mate the restitution coefficient of our system. We obtain 
e = 0.95 ± 0.03, which compares favorably with values 
known from the literature 

Intermediate regime: As more liquid is added, the 
width of the hysteresis loop becomes larger, such that the 
kinetic energy of the impact may sometimes not be suffi- 
cient to break a bridge. It thus appears likely that small 
clusters will form which have a certain lifetime. Since 
these are heavier than the single beads, their velocity af- 
ter an impact is less affected by the liquid bridges. Fur- 
thermore, the interior of a (larger) cluster may act as a 
sink for the liquid due to capillarity, resulting in a some- 
what reduced wetness of the beads outside the cluster. 
Hence the 'effective' restitution coefficient, averaged over 
all entities in the system, increases again. As a conse- 
quence, the BNE is now reduced, and the phases become 
better miscible, as we observe. 

This model is in accordance with the observation that 



the overall appearance of the material is not expected 
to change noticeably in this regime: as more liquid is 
added, the average cluster size is just increased. The 
clusters take up the excess liquid, due to the increased av- 
erage coordination number of the bridge network within 
a cluster. In other words, excess liquid is stored in liq- 
uid bridges within clusters, which are not active in the 
dissipation process. 

We thus have elucidated the interplay between segre- 
gation (BNE) and a phase transition to viscoplasticity, 
and have given a robust criterion for the latter. The 
phase boundary can be understood from simple energy 
considerations. 
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